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PREFACE 


The purpose of this research was to investigate the utility of 
Rossby's wave equation in the westerlies at the 600 millibar level by 
statistioal methods. Also a statistical analysis was made to investigate 
just how wave length, tonal index and latitude contribute to the move¬ 
ment of trouts in the westerlies* 

This work was conducted at the Chited States Haval Postgraduate 
School, Monterey, California, during the period December 1948 to May 
1949* It was done to moot partial requirements for the degree of Master 
of Soi anoe in Aerology* The idea for the investigation was suggested 
by Professor F* L* Martin of the Staff of the Department of Aerologioal 
Engineering of the United States Naval Postgraduate School* His 
assistance and advice in the preparation of this thesis is gratefully 
acknowledged* 
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SYMBOLS AND ABBREVIATIONS 

v Eastward speed of trough at 500 millibars as computed by Rosaby’s 
ware equation* 

I - Actual eastward speed of trough at 500 millibars* 

* Mean west to east wind component* 

- Height difference for 10° of latitude, 5° on each side of mean 
latitude. 

- Northward rate of change of Coriolis par octet er. 

Actual wave length, measured from the trough involved westward to 
the next major trough, 

- Stationary wave length* 

- Latitude* 

- Mean radius of earth. 

- Angular velocity of earth. 

- Standard deviation* 

- Correlation coefficient. 

■ i '. • 

- Standard error of estimate* 



CHAPTER I 


INTRODUCTION 


The problem of forecasting the movement of long waves in the westerlies 
is very important in Meteorology. Although several approaches and techniques 
have been developed in an attempt to forecast the movement of upper air 
troughs, no perfect technique has been developed to accomplish this aim* 

As yet, no theoretical or statistical investigation has shown that the move- 
ment of these upper air troughs oan be predicted with any great degree of 
accuracy* 

The purpose of this paper is to sake an objective Investigation of 
the forecasting of the movement of major troughs in the westerlies at the 

•V' • » 

600 millibar level* The movement of these major troughs will be investi¬ 
gated in the zone 30 to 60 degrees north latitude, mainly over the continent 
of North Amerioa and adjaoent oceans* The problem is complex and it will 
be approached in several ways, mainly statistical* 

It will be of interest to consider previous approaches developed in 
an attempt to solve this problem. First, there is the purely kinematio 
approach of S* Potterssen (11)* We will, however, concern ourselves pri¬ 
marily with tho dynamical and statistical approach* 

The main dynamical approach is that of Rossby (12) who developed the 
equation 
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This equation gives the west-east velocity, , of sinusoidal waves in 
the westerlies, where U is the speed of the westerlies (assumed eonstant), 

(1) 




J@ is the northward rate of change of the Coriolis paroaeter and L is the 
ware length* In addition to the above-mentioned assumptions, Rossby also 
assumed a homogeneous incompressible atmosphere, zero horizontal divergence, 
and a scale of motion large enough so that tho variation of the Coriolis 
parameter with latitude is of importance in determing wave dimensions* 

The theory Was extended by Hauradtz (6) to take into account the curva¬ 
ture of the earth and the finite lateral extent of the wave. Holaboe (8) 
obtained similar results for the more general barotropio atmosphere. 

In 1944, J. Bjerknes and Holaboe (1) presented a theory of wave motion 
in a baroclinio atmosphere. The field of horizontal divergence is evaluated 
from the pressure pattern by means of gradient-wind relationships, and on 
this basis the following relation is established* If denotes the 

speed of the westerly current at any height <4£c an increasing function 
of wave length and o the wave velocity, then 

<2) C « sCL(h) - -66 c 

where h ie the height at which the mass divergence in the horizontal velocity 
field is zero* 

Charney (3) considered the case of a baroclinio atmospheric model having 
a mean zonal wind In agreement with that in middle latitudes for the month 
of February* This model, however, has no horizontal wind shear* It is shown 
that Sossby’ s simple wave formula for the speed of a barotropic wave will 
apply to the baroclinio wave if the constant value of the zonal wind in the 
formula is the mean zonal wind averaged with respect to pressure from the 
top to the bottom of the baroollnie atmosphere* For the February model 
already mentioned, Charney finds that this mean zonal wind Is that in the 
vicinity of the 600 mb* level* It has bean pointed out by Holaboe that the 
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formula for the barocllnio ware speed is strictly true only at the level of 
non-divergence in the atmosphere, which level, therefore, is in the vicinity 
of 600 mb* 

Due to the fact that Rossby’s original formula (1) was derived on the 
basis of sinusoidal wives in the westerlies, it suggests the importance of 
applying it to major atmospheric waves rather than minor waves which disappear 
at relatively low levels in the atmosphere. Hernias and Clapp (10) used 
Rossby’a formula for calculating the movement of five day mean troughs at 
45 degrees latitude during the winter. The correlation between this calcu¬ 
lated movement and the actual full week trough movement was .66 when cases 
of new trougjh development were eliminated. In contrast to this, the appli¬ 
cation by Brooks (2) and Yehle (14) to daily trough displacements yielded 
very poor results, but they apparently neglected to differentiate between 
major and minor waves. However, Craig subsequently obtained a much better 
correlation by using major waves. This was done by selecting as the wave 
length the distance between the trough concerned and the next one to the 
west having approximately the same depth. Possibly, this improved result 
was due to the fact that seleotion of the most pronounced troughs is com¬ 
parable to smoothing pressure patterns by constructing five day means. 

Finally, Cressman (4) made an important contribution to the statistical 
study of long wave movements in the westerlies* He expressed Rossby’s 
formula in the form 

w C 5 (Is - L*) 

where Ls is the wave length of a hypothetical stationary wave defined by 

/Of *~ 

the mean tonal wind, B v —• Then the trough velocity can be determined 
by either Ls or B. Cressman constructed for each of 120 individual cases 
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a meridional profile of the tonal wind at 600 mb., the level of non- 
divergence. This enabled him to select the proper value of 0 at each 
latitude. The latitude to be used was arbitrarily selected as the 
latitude of the maximum west wind on the profile. The meridional profile 
of the zonal wind at 600 mb, was computed with the aid of the hydrostatic 
equation from meridional profiles of the zonal wind and the temperature 
gradient at 700 sib. 

By testing the wave formula and using the above technique, Cressman 
found an absolute error of only two degrees of longitude per day* However* 
he selected for study those special synoptic situations when the major waves 
were qxmsistati onary or were moving with a constant speed for a period of 
several days. 

Cressman*s method depends upon a laborious extrapolation process te 
the 600 sto* level and thus is not a practical tool. In order to reduce 
some of this labor* he attempts to determine L g from the 700 mb. chart. 

This is done by computing constant vortieity trajectories eentered on the 
700 sib. jet stream. He found that a value of lag determined by using a 
700 mb. mean zonal wind tJ centered on the jet stream and averaged over 4° 
latitude on either aide of the jet, gives good agreement with the value of 
L s determined from the 600 mb. zonal wind profile. By this device, the 
problem of determining L g is brought back by objective methods to the 
700 mb. level. 

The results obtained by Cressman using (3) together with the value of 
Lg obtained from the 700 mb. chart, does not give as favorable results as 
the method based on the 600 mb. chart. In fact, he states, "the data showed 
no simple relation between the speed of the wave and the difference between 
stationary and actual wave lengths," Thus, while the 600 nb, method of 
analysis gives fair results, it is too laborious, while methods based on 
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tho use of the 700 mb. chart arc not accurate enough for normal verifi¬ 
cation. 

Even though the level of non-divergence is between 610 and 670 mb., 
according to Charney, it was decided to test Rossby*s formula at the 600 mb. 
level. While it is expected that the value of TJderived from this chart 
will b© too large, it is believed that it may be consistently too large, 
thus giving a good correlation to the desired correct value of 0 from the 
600 mb. chart. 

Tho data used in this investigation were obtained from northern 
Hemisphere Historical Weather Haps (7) which are drawn for the surface and 
the 600 mb. level. These maps, published quite recently, are considered 
the best Northern Hemisphere charts available at this time. The charts 
for the months of February, March and April, 1946, were used for this 
particular paper. 

Taro statistical methods of attack are pursued, namely* 

1. A strict teat of the Bossby formula using 600 mb* data by comparing 
C* with actual trough speed C. 

2. A more general statistical analysis of the problem based on determin¬ 
ing the 24-hour trough velocity C as a linear regression on the 
variates, wave length L, the mean latitude^?, and the height dif¬ 
ference n € . 

Tho solution of problems 1 and 2 necessitates objective methods of 
determining the variates. These objective methods are described in the 
following chapter* 


(5) 



CHAPTER II 
METHOD OP AHALY3IS 

This chapter will describe in detail the method of analysis need in 
this investigation* It will consist mainly of the procedure used la 
obtaining the variates from the synoptic charts at the 500 mb* level* 

This procedure will also be the basis for using graphs based on this 
method of analysis* 

Mnor waves have their greatest amplitude in the low levels and only 
very small amplitudes at the height of tho tropopause* and in general move 
very rapidly* However* major waves may be almost obsoure at sea level* but 
the amplitude of these waves increases with height and reach their greatest 
strength at the tropopause, in this investigation, only major waves will 
be used* a major wave being defined as one which has a well-defined ridge 
lying between two adjacent troughs* 

The wave length of these major waves is almost always so large* that 
if one trough is over the dense network of observations of Forth America* 
the other troughs are over the Pacific or Atlantic Ooean* Due to this fact* 
it is usually Impossible to locate an upper air trough closer than to 
within five degrees of longitude, if this is true* it is reasonable to 
assume that the expected error in determining the actual wave length would 
be about seven percent if m average wave length of 70 degrees longitude is 
assumed* 

Troughs were considered to be at the point of maximum curvature of the 
contour lines* but where the curvature was constant* the southernmost point 
of the oontour line was used* 

For utility* the wave formula may be given in a more convenient notation* 
following the convention that velocities are to be measured in degrees of 
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longitude over which the wind or pertubatioa moves in & period of twenty- 
four hours and the wave length of the pertubation is measured in degrees 
of longitude rather than linear measure. Thus equation (1) becomes 




where <y is the latitude* 

The objective methods of selecting variates are as follows* 

1* Mean latitude of the wave* A contour line was selected which 

crossed the trough involved and the next trough to the west of it* 
This contour line was imbedded in the mean westerly current of the 
trough involved. Such a contour line is shown in Figure 1 as T5f' 
The intersection of this contour line with the.trough lines gave 
latitudes <P, t and * The average of ^ , and ^ gave the mean 
latitude* In this investigation, this mean latitude was rounded 
off to the nearest five degrees, probably without causing any 
error greater than observational errors* 

2* The wave length was measured as the distance in degrees longitude 
from the trough under investigation westward to the next trough at 
the mean latitude* 

3* The mean zonal wind was measured in a five degree tone both sides 
of the mean latitude and bounded by the trough involved and the 
next trough to the west* 

4* Modification of Rossby*s formula to constant pressure analysis* 
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A Z is the average height difference per ten degrees of latitude 
centered on the mean latitude. This modification is based upon 
the assumption that the sonal "wind is the geostrophic wind* 

In seleoting synoptic situations, there were only two limitations* 

First, there must be a broadband of westerlies which did not contain any 
closed centers, since the east winds associated with the closed centers would 
tend to cancel out an equal amount of west winds in averaging out the height 
differences oo that the final sonal wind profile would not give an accurate 
representation of the strength of the westerlies» Secondly, only waves Which 
approached sinusoidal curves, end did not have exceptionally large or small 
amplitudes were used* 

In all, 88 oases satisfying these criteria were selected during the 
months of February, Maroh and April, 1946, with a mean latitude range of 
50 to 50 degrees north latitude* 





CHAPTER III 

RESULTS AND INTERPRETATION 

All values of speed shown in these results will be in units of degrees 
of longitude per day. The results of a test of Rossby’s formula follows 

c * jo.y. ± ifij 

i 

% = //.?■£ 

The linear regression of C on C is given by 

( 6 ) C - .31C* * S. 78 

with a correlation coefficient between C and C A of .49* The standard error 
of estimate in estimating C from was 3.6P longitude per day. 

From the mean values shown above, it is seen that the wave formula 
tends on the average to overestimate the speed of the troughs by a degree 
and a half. This is not surprising, sinoe strictly speaking, the formula is 
valid only in the vicinity of the non-divergence level, approximately 600 mbs. 
However, from the regression equation, C and C R ere equal when they have a 
value of 9.46 degrees longitude per day* Therefore, the conclusion already 
drawn from the arithmetic means alone is only true when C* is greater than 
9.46, which is the case in the majority of observations. When C* is less 
than 9.46, C is greater than the calculated C A • 

Now we will correlate the variates C on d 2 , L and <P • The units are 
as follows t C in degrees longitude per day, d a in hundreds of feet, L in 
degrees longitude, and <P in radians. The simple correlations are shown in 

* Standard deviation was used here. See Kenney (9), p. 114* 
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Table II• The results of the multi-variate analysis followt < 

Multiple correlation coefficient 

(c-**,L;<P) * .63 

Standard error of estimate 

=3./3' L oh yt^c/e/Day 

linear regression equation 

<n c = JGAZ -.ni +!T16f A/. 08 

In a test of regression equation (7) the following distribution of 
computed minus observed displacements are shown here* The units are degrees 
longitude per day. 

< - 6° -6° to -2° -2° to +2° 42° to *6° >6° 

3 12 33 16 2 

It is seen the regression equation (7) gives a normal distribution with 
the errors distributed almost equally on each side of the oentral block. In 
62.3^ of the oases, the difference is less than 2 degrees longitude end 94% 
of the cases show a difference of less than 6 degrees* 

The average absolute errors as given by regression equations (6) and (7) 
are 2.76 and 2.46 degrees longitude por day, respectively. 

In interpreting the results, it iB seen that the simple correlation 
coefficients agree qualitatively with Hossby’s wave formula (1), The corre¬ 
lation coefficient of-f>%34 between the observed speed of the trough and the 
height difference indicates that the speed of the trough varies directly with 
the zonal index. Also the correlation coefficient of -.30 between the speed 
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TABLE I 




Means of Variates 


c 

10,41 

A & 

6.38 

L 

45.42 

<P 

•70 

Standard Deviation 


4.1 

1.4 

11.2 

.10 


c 

— 

•34 

-.30 

.53 

Simple . _ 

✓ 

41 It 

.84 

. — 

•05 

•68 

Correlation . 


L 

-.SO 

u> 

o 

• 

— 

.06 

Coefficients 


•63 

• 58 

.06 
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of the trough end the wave length shows that the speed of the trough de¬ 
creases as wave length increases. 

Probably of greater significance is the fact that a correlation coeffi¬ 
cient of +.53 was found between the speed of the trough and the latitude* 
this effect of latitude agrees quite well with Eosaby * a wave formula (l) 
where the value of $ decreases by about 26$ from SO to SO degrees latitude* 
Although the latitude has an effect on both terms of equation (5), its 
effect on the last term of the equation is more predominant, which agrees 
well with the results* - - 

The simple correlation coefficient of f .68 between the latitude and 
the height difference is due to the closed highs at low latitudes as shown 
by the synoptic charts* These closed highs oause a low zonal index which 
is shown in the results* 

Although the simple oorrel&tion coefficients agree quite well with 
Hossby’s wave formula (1) qualitatively, they do not agree quantitatively* 
This is shown by the fact that the correlation coefficient between C and 
Cr is only *49, with a standard error of estimate for C of 3.6° longitude 
per day* 

By developing the linear regression equation (7), the correlation 
coefficient was raised to *63. While this is considerably better than *49, 
it does not lower the standard error of estimate appreciably* However, 
equation (7) does help to show quantitatively the effect of height difference 
(A2-), wave length (l»), and latitude {(f) on the speed of the trough (0)* 

For an average case, the latitude effect is about 12 times as important as 
the height difference term and about twice as important as the wave length 
term* This would indicate that at a specific latitude, the movement of the 
trough would depend primarily on the wave length and that the zonal index 
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plays a very small part. Ia fact, the tonal index term would have a maximu m 
value of 1.6° longitude per day for the average case. With present methods 
of obtaining meteorological information, the movement of a trough at 600 sib* 
could not be measured to 1 or 1.6 degrees oif longitude, so the zonal index 
term may be considered negligible as compared to the wave length and lati¬ 
tude terms. These results agree quite well with those of Wandas and Clapp- 
( 10 ), 

As the development of a tool for the forecasting of movement of troughs 
at the 500 mb. was one of the primary objectives of this paper, a graph ia 
plotted which gives trough velocity when the mean latitude, height difference 
and wave length are known. However, it is necessary to use the technique 
of analysis as described in Chapter II. The main conclusions are stated 
in the following* 

SOKMART 

1. Rossby's wave equation (1) is not applicable at the 500 mb. level. 

2. Equation (6), when used to estimate trough speed, gives a standard 

error of estimate of 5.6° longitude per day with an average absolute 
error of 2.76° longitude per day* 

3. Equation (7), when used to estimate trough speed, gives a standard 

error of estimate of 3,13° longitude per day with an average abso¬ 

lute error of 2.46° longitude per day. 

4. The zonal index plays a relatively small part in the movement of 
troughs at the 600 mb* level. This is further substantiated by an 
analysis of the synoptio 600 mb. ohart. In every case where there 
was a closed high behind the trough, the observed movement was 
greater than the calculated movement. This indloates that the 
trough moved even though there was a very weak tonal flow from west 
to east. 

(14) 


The fact that the movement of the trough as calculated by Rossby’a 
wave equation (1) overestimates the actual trough movement indi¬ 
cates that the level of non-divergence is below 600 mb* 

Even after considering. the errors introduced due to present observa¬ 
tional techniques, it is apparent -that there must be another variate 
which contributes to the speed of the trough. It is quite possible 
that this variate is wave amplitude* 
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APPENDIX I 
VARIATION CP& WITH LATITUDE 


90® 

76° 

60° 

50 ® 

45® 

30 ® 

15® 

0 ® 


0.0 

0.593 

1.145 

1.461 

1.619 

1.983 

2.212 

2.290 




